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ABSTRACT 

Dust grains grow in interstellar clouds by accretion and coagulation. In this paper, we focus on 
these two grain growth processes and numerically investigate how they interplay to increase 
the grain radii. We show that accretion efficiently depletes grains with radii a < 0.001 |J.m 
on a time-scale of < 10 Myr in solar-metallicity molecular clouds. Coagulation also occurs 
on a similar time-scale, but accretion is more efficient in producing a large bump in the grain 
size distribution. Coagulation further pushes the grains to larger sizes after a major part of 
the gas phase metals are used up. Similar grain sizes are achieved by coagulation regardless 
of whether accretion takes place or not; in this sense, accretion and coagulation modify the 
grain size distribution independently. The increase of the total dust mass in a cloud is also 
investigated. We show that coagulation slightly 'suppresses' dust mass growth by accretion 
but that this effect is slight enough to be neglected in considering the grain mass budget in 
galaxies. Finally we examine how accretion and coagulation affect the extinction curve: The 
ultraviolet slope and the carbon bump are enhanced by accretion, while they are flattened by 
coagulation. 

Key words: dust, extinction — galaxies: evolution — galaxies: ISM — ISM: clouds — ISM: 
evolution — turbulence 



1 INTRODUCTION 

Dust enrichment in galaxies is one of the most important top- 
ics in galaxy evolution. Dust grains actually modify the spec- 
tral energy distribution of galaxies b y reprocessing stellar radia- 
tion i nto far-infrared wavelengths (e.g. lDesert. Boulanger. & Pugell 
Il99d) . Grain surfaces also regulate interstellar chemistry; espe- 
cially, formation of molecular hydrogen predominantly occurs on 
dust grains if the interstellar medium (ISM) is enriched with dust 
(e.g.'Cazaux & Tielens These effects of dust grains point to 

the im portance of clarifying the dust enrichment in galaxies (see 
lYamasawa et al.i201 ll for a recent modeling). 

Dust enrichment is gov erned by vari ous processes depend- 
ing on age, metallicity, etc. jDwekl Il998h . Dust grains are sup- 
plied by stellar sources such as supernovae (SNe) and asymptotic 
giant branch (AGB) stars (e.g. Kozasa et al. 2009; Valiante et al. 
l2009h . Dust is also destroyed by SN shocks in the ISM. In the 
Milky Way, the ti me-scale of dust destruction by SN shocks 
is a few x 10^ yr jjones. Tielens. & HollenbachI 1 19961 : but see 
I Jones & Nuthll201 ll) . while that of dust supply from stellar sources 
is longer than 1 Gyr jMcK ee 1989). Therefore, to explain a sig- 
nificant amount of dust in the ISM, it has been argued that dust 
;rains grow i n the ISM by the accretion of metal s onto grains 
Inoud 120031 : iDraind I2OO9I: IZhukovska elal] I2OO8I : iPipino et all 



E-mail: hirashita@asiaa.sinica.edu. tw 



I2OIII : IValiante et al.ll2oTll: lAsano et alj|2012l) . Dust grains grow 
efficiently in molecular clouds, where th e typical number density 
of hydrogen molecules is ~ 10 cm" fairashitallioooh . Larger 
de pleti on of metal elemen ts in cold clouds than in warm medium 
JSavage & Sembaclj|l996h may indicate grain growth in clouds. 

In dense environments, dust grains grow not only by accre- 
tion but also by coagulation. Indeed, a deficit of very small grains 
contributing to the 60 vim emission is observed around a typical 
density in molecular clouds > 10^ cm~^, and is interpreted to be a 
consequence of coagulation ( Stepnik et alj|2003l) . Thus, we should 
treat accretion and coagulation at the same time, since these two 
processes could compete or collaborate with each other. Moreover, 
accretio n and coagulation may have differ ent impacts on the ex- 
tinction dCardelli. Clayton. & Mathisll 19891) : Coagulation shifts the 
grain sizes to larger ranges and flattens the extinction curve, while 
accretion also increases the extinction itself. Which of these two 
effects dominates can also be clarified by solving accretion and co- 
agulation simultaneously. Thus, we examine in this paper how these 
two major processes of grain growth - accretion and coagulation - 
interplay. We also investigate the impact of grain growth in clouds 
on the extinction curve. 

When the dense clouds are dispersed after their lifetimes, the 
dust grains grown through accretion and coagulation are injected 
into the diffuse medium, where dust destruction by SN shocks oc- 
curs CMcKee. 1989) . Hirashita & Yan (2009) show that grain shat- 
tering is also efficient in the diffuse ISM. Thus, our calculation re- 
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suits in this paper, which focuses on grain growth in clouds, can 
be used as 'inputs' for subsequent grain destruction and shatter- 
ing. SN shocks do not penetrate efficiently into dense environ- 
ments and dust destruction by SN shocks can be neglected in dense 
clouds jMcKee 1989). Other possible dust processing mechanisms 
in dense clouds such as destruction by protostellar jets and photo- 
destruction by radiation from stars are also neglected, since the ef- 
ficiencies of these processes are not clear yet. Thus, we do not treat 
these destructive processes and concentrate on the grain growth 
mechanisms (i.e. accretion and coagulation) that work in the dense 
ISM. 

This paper is organized as follows. We explain the formula- 
tion in Section |2l and describe some basic results on the evolution 
of grain size distribution through grain growth (accretion and coag- 
ulation) in individual clouds in Section[3l Based on the calculation 
results, we discuss effects on the extinction curve and implication 
for the galaxy evolution in Section |4] Finally, Section |5] gives the 
conclusion. 



2 FORMULATION 

In this section, we formulate the evolution of grain size distribution 
by the accretion of metals and the growth by coagulation in an in- 
terstellar cloudQ These two processes are simply called accretion 
and coagulation in this paper. We particularly focus on the interplay 
between accretion and coagulation. 

Throughout this paper, we call the elements composing grains 
'metals'. We only treat grains refractory enough to survive after 
the dispersal of the cloud, and do not consider volatile grains such 
as water ice. We also assume that the grains are spherical with a 
constant material density s, so that the grain mass m and the grain 
radius a are related as 
4 3 

m = -7ra s. (1) 
3 

We define the grain size distribution such that n{a, t) da is 
the number density of grains whose radii are between a and a + da 
at time t. For simplicity, we assume that the gas density is con- 
stant and that the evolution of grain size distribution occurs only 
through accretion and coagulation. The dust mass density, pd it) is 
estimated as 



Pd(t) 



^ s n{a, t) da. 



(2) 



We adopt silicate and graphite as dominant grain species (e.g. 
iDraine & Leell984h . and treat these two species separately to avoid 
the complexity arising from compound species. 

Accretion and coagulation are separately described in the fol- 
lowing subsections, but we solve these two processes simultane- 
ously in the calculation. The details of the numerical schemes are 
explained in Appendix IaI 

2.1 Accretion 

Because t he knowledge about chemical properties of accretion is 
still poor jjones & Nuthll20T l|), we simplify the picture by assum- 
ing that grain growth is regulated by the sticking of the key species 

^ Although we mainly consider a molecular cloud for a 'cloud', the for- 
mulation is not specific to molecular clouds but applicable to any clouds 
including cold neutral clouds. Therefore, we simply call the place hosting 
grain growth 'cloud' . 



denoted as X (X is Si a nd C for silicate and graphite, respectively; 
lHirashita&Kuoll201ll hereafter HKll). The mass fraction of the 
key species in dust is denoted as fx- fx ~ 0.166 for silicate (i.e. a 
fraction of 0.166 of silicate is composed of Si), while /x = 1 for 
graphite (i.e. graphite is composed of only C) (Table [T|l. We follow 
the formulation in our previous paper (HKl 1 ; see also Evansll994h , 
but modified for the purpose of numerical calculations. 

Since the grain number is conserved in accretion, the follow- 
ing continuity equation in terms of n{a, t) holds: 



dn(a, t) dr. , 



0, 



(3) 



dt 

where a = da/dt is the growth rate of the grain radius, which is 
given by the following form (HKl 1): 



= C(t) a/r(a) 



(4) 



Here the growth time-scale as a function of grain radius, r(a), is 
estimated as 



r(a) 



nx.tot mxS ( kTr 



fxs \2nmxJ 



1/2 ■ 



(5) 



where nx,tot is the number density of element X in both gas and 
dust phases, mx is the atomic mass of X, S is the sticking prob- 
ability for accretion, fx is the mass fraction of the key species in 
the dust (see above), k is the Boltzmann constant, and Tgas is the 
gas temperature. In equation we also introduce the fraction of 
element X in gas phase: 



^{t) = nx(t)/nx,tot, 



(6) 



where nx{t) is the number density of element X in gas phase as 
a function of time. Note that d is independent of a. The gas-phase 
metals decrease by accretion as 



dnx(f) 
dt 



47ra'nx (t) f -^^^^ S n{a, t) da 
\ z-rnnx / 



(7) 



Since we are often interested in the grain mass, it will be con- 
venient to consider the grain size distribution per unit grain mass 
rather than per unit grain radius. Thus, we define n(m, t) as the 
number density of grains with mass between m and m + dm. The 
two functions, h and n, are related by n(m, t) dm = n{a, t) da; 
that is, n — n/ (47ra'^s) by using equation (TJ. Then, the time evo- 
lution of h is obtained from equation ^ as 



dimn) , d{mn) 1 . . 



(8) 



where we have used equation l[T), m = dm/dt = 47ra^sd, 
d/dm = [l/(47ra2s)]a/9a, and dh/da = 0. 
It is convenient to define a as 



(j(m, t) = mn{m, t). 

Then, equation ^ is reduced to 

do . da 1 . 



(9) 



(10) 



where /i = In m and fi = d/i/dt . Noting that /i = 3a /a, we obtain 
from equation Q 



(11) 



r(m) ' 

where r is now expressed as a function of m instead of a. The 
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evolution of ^ is calculated by (equations |5] (6) [T] and[9} 



At 



mxnx.tot 



■ dm. 



(12) 



10 -rim) 

We solve equations JlOt . Jilt , and J12b . The discretized form of 
equation dlOt is shown in Appendix lAll The initial conditions are 
described in Section|23] 



2.2 Coagulation 

For the evolution of grain size distribution by coag ulation, we apply 
the fo rmulation developed in our previous work jHirashita & Yan" 
120091) . Here we adopt an analytic formula derived by Ormel et al. 
1 1200 9) for the turbulent velocities. Below we briefly overview the 
method, focusing on the treatment of the grain velocities. 

We consider thermal (Brownian) motion and turbulent mo- 
tion as a function of grain mass (or radius, which is re- 
lated to the mass by equation [TJ. The turb ulent motion be- 
comes important particularly for large grains tossenkopj Il993l : 
IWeidenschilling & Ruzmaikinalll994| - lormel et al...2009.) . The ve- 
locity as a function of grain mass m is given by a combination of 
thermal (Brownian) and turbulent velocities (uth and uturb, respec- 
tively) as 

v(mf = vti,{mf + ■uturb(m)^. (13) 

The t hermal velocity is given by v^y^ = 8A;rgas/(7rm) jSpitzeil 
Il978t) . which is numerically evaluated as 



0.529 



1/2 



0.1 ^m 



-3/2 



3 g cm~ 



-1/2 



(14) 



The velocity driven by turbulence is given by 



= 8.8xl0^(g|) 



1/4/ „ 



0.1 ^m 



VIO^ 



^mol 



-1/4 



3 g cm 



1/2 



cms ^, (15) 



where rimoi is the number density of the molecular gas, which is 
related to the number density of hyd rogen nuclei, nn, as rimoi = 
nn/l.T for the cosmic abundance dOrmel et alj|200^jfl In fact, 
equation dlSt represents the relative velocity of equal-sized grains. 
We adopt this as a typical velocity for the following reasons: (i) 
a simple form for the velocity as a function of grain size fits our 
approximate treatment of the relative velocity (see below); (ii) the 
relative velocity is determined by the velocity of the larger grain, 
which is coupled with the larger-scale turbulent motion, so the rel- 
ative velocity is always of the order of i'turb(t[i), where a\ is the 

^ 'T^mol = l-7nH in lOrmel et alj i2009l) should be n^oi = nn/l.T 
(note that n,jjol is denoted as n in lOrmel et alj|2009h . We can also eas- 
ily chec k the validity of t he "intemiediate regime" by using the esti- 
mates in lOrmeletai] j2009h : with Tgas = 10 K and nn = 10^ cm-^, 
St ~ 1.3 X 10-2 (a/0.1 nm) and Rc^^/^ ^ 4 g ^ io-4, where St is 
the Stokes number and Re is the Reynolds number. Thus, the condition for 
the intermediate regime. Re" ^''^ < St < 1, is satisfied in the size range 
where the turbulent motion is dominant over the thermal motion. Yan et al. 

adopted larger size and velocity for the largest eddies, so Re is larger 
than assumed above. In such a case, the intermediate regime is completely 
valid for all grain sizes treated in this paper 



radius of the larger grain; (iii) the uncertainty caused by this rough 
treatment of relative velo cities, compared with th e analytic solution 
given by equation (28) o f lOrmel & Cuzzil j2007h . is within a factor 
of 1.3. 

Grain velocities are dominated by thermal motion for small 
grains (a < 0.002 iJ.m), while turbulence drives the motion of large 
grains efficiently. If the grain size is too large, the grain velocity 
becomes larger than the coagulation thresho ld, which is calculated 
by the same way as iHirashita & Yaiil ( l2009l) . Yet, the grain veloc- 
ities ( ^ 1 km s"^) a re too small for shattering or erosion to oc- 
cur ( Jones etalj|l996h . For a grain colliding with a ^ 0.01 iJ.m 
grains, coagulation occurs if the grain has a size smaller than 
a few X 0.01 \im. Each time-step is divided into four equal small 
steps, and we apply Vki = Vk + ve, \vk — vi\, Vk, and vi {vk and 
Vi are the grain velocities evaluated at a = and ag, respectively, 
and Vki is the relative velocity in the collision between two grains in 
bins k and £; see Appendix IA2l for the discrete formulat ion) in each 
step to consider a variety of r elative velocity directions ( I Jones et al] 
1 1 994l; IHirashita & Yanll2009l) . 



2.3 Initial conditions 

The initial grain size distribution is assumed to be described by 
a power-law function with power index ~r and upper and lower 
bounds for the grain radii Omin and flmax, respectively: 



n{a, 0) 



(4-r)pd(0) 



4 / 



4 — r 
max 



(16) 



for Qmin < a ^ Hmax. If a < Qmin Or a > Qmax, 0) = 0. 

The dust mass density at t = 0, pd(0) satisfies equation l|2ll, and is 
given later by equation dlSt . The dust mass density is related to the 
initial condition for ^ as pd(0) — (mx//x)[l — ^(O)]nx,tot, and 
the total number density of element X both in gas and dust phases 
is written as 



7^x,tot 



^0 



(17) 



where Z is the metallicity, and (X/H)q is the solar abundance (the 
ratio of the number of X nuclei to that of hydrogen nuclei at the 
solar metallicity), and tih is the number density of hydrogen nuclei. 
Therefore, pd(0) is related to the initial condition for ^ as 



Pd(0) = 



(18) 



We apply Z — Zq (note that the time-scales of both accretion 
and coagulation are simply scaled as Z~^). We assume ^(0) — 
0.3, which roughly match es the depletion in the diffuse medium 
dSavage & SembachI 19961) . Although there is uncertainty in the de- 
pletion because of the assumed elemental abundance pattern (usu- 
ally the solar abundance pattern is assumed), the following discus- 
sions on the evolution of grain size distribution are not altered as 
long as we adopt a reasonable value such as ^(0) = 0.1-0.7. 

iMathis. Rumpl. & Nordsiec3 ( Il977h show that the extinction 
curve in the Milky Way can be fitted with a power-law grain size 
distribution with r = 3.5. Thus, we assume that r — 3.5. The effect 
of r on accretion has already been investigated by HKII. Briefly, 
for larger r, small grains occupy a larger fraction of total dust sur- 
face and the total grain surface becomes larger; as a consequence, 
the grain growth occurs more rapidly for larger r. Since the largest 
grains are not susceptible to accretion and coagulatio n as we show 
later, we fix the maximum size as ttmax ~ 0.25 ^.m jMathis et al.l 



© 2012 RAS, MNRAS OQO.fTlflOl 



4 Hirashita 



The lower bound of the grain size is poorly determined from 
the extinction curve (Wei ngartner & Draine.2001i) : thus, we exam- 
ine ttmin = 0.3 and 1 nm. 

2.4 Selection of quantities 

As mentioned at the beginning of this section, we consider sili- 
cate and graphite separately. The quantities adopted in this paper 
are summarized in Table[T] and are based on HKll. Since we are 
interested in the interplay between accretion and coagulation, we 
fix parameters which do not affect the relation between those two 
processes. 

By using equation dlVt . r(a) in equation (O can be estimated 

as 



Table 2. Models. 



r = 1.61 X 10" 



0.1 [im 



riH 



103 cm-3 



\wk) Vo.sy 



(19) 



for silicate, and 



= 0.993 X 10'' 



0.1 |J.m 

1/2 

0.3 



103 : 



(20) 

Z = Zr 



,10K 

for graphite. As mentioned in Section [23] we adopt zj = zjq 
riH = 10^ cm~'^ for the typical values derive d from observa - 
tional properti es of Galactic molec ular clouds ( Hi rashital |2000|) . 
Tg^, = 10 K dWilson. Walker, & Thornlev 1 997), and S = 0.3 
jLeitch-Devlin & Williamsll 19851 : lOr'assi et al.ll201 ih . Both accre- 
tion and coagulation have the same dependence of time-scale on 
metallicity and density as oc (nn^) ^^ 

For t, we in vestigate a probable range for t he lifetime of 
molecular clouds. iLada, Lombardi, & AlvesI ( I2OI0I) mention that 
molecular clouds survive after star formation activities lasting 
2 Myr. The comparison with the age of stellar clusters associated 
with mo lecular clouds indicates that the lifetime of clouds is ~ 
10 M yr iLeisawitz, Bash, & T haddeusI 1 19891: iFukui & Kawamur j 
l20ld) . iKoda et al. ( 2009 ) argue that molecular clouds can be sus- 
tained over the circular time-scale in a spiral galaxy (~ 100 Myr). 
Therefore, we examine t ~ a few~100 Myr as a probable range for 
the cloud lifetime. 



3 RESULTS 

3.1 Evolution of grain size distribution by accretion and 
coagulation 

We examine various cases, focusing on the interplay between ac- 
cretion and coagulation. We fix r and Omax since they do not affect 
the relative role between accretion and coagulation (Section [3.31 ). 
Other parameters, especially amin and turbulent velocity, are po- 
tentially important in this paper. For convenience, we name each 
set of parameters Model A-D as shown in Table|2] 

In Fig.[T] we present the results for Model A alt = 10 Myr. 
To show the mass distribution per logarithmic size, we multiply 
to n. Since the growth rate of grain radius, d, is independent 
of a, the impact of grain growth is significant at small grain sizes. 
Moreover, gas-phase metals accrete selectively onto small grains 
because the grain surface is dominated by small grains (see also 
IWeingartner & Draine|[T999l) . The results for silicate and graphite 



Model 


r 






turbulence 






[nm] 


[^m] 




A 


3.5 


0.3 


0.25 


on 


B 


3.5 


1 


0.25 


on 


C 


3.5 


0.3 


0.25 


off 


D 


3.5 


1 


0.25 


off 



are similar because the amount of available gas-phase metals is 
similar. 

In order to clarify the contributions from accretion and coag- 
ulation. Fig. [T] also shows the cases where either only accretion or 
only coagulation is taken into account in the calculation. In Fig.[T] 
we observe that both accretion and coagulation deplete grains at 
a < 0.001 \im. Accretion plays a significant role in increasing the 
grain mass around a ~ 0.001 [im, and coagulation further pushes 
the peak to a ~ 0.002 |i.m at 10 Myr. 

Fig. [T] also presents the analytical solution for accretion cal- 
culated by the method in HKll. Our numerical results with only 
accretion match the analytical solution very well except for the 
peak of the size distribution around a ~ 0.001 \im. The sharp peak 
is smoothed because of the numerical diffusion. Reproducing the 
sharp peak is not important, however, since it is smoothed out by 
coagulation in any case. 

Next, we show the results for Model A at various t in Fig. [2] 
Fig. [2^ presents the grain growth by accretion and coagulation. 
We observe continuous grain growth. The spiky features around 
a ~ 0.01-0.06 [im are produced by the accumulation of grains at 
sizes where the grain velocity reaches the coagulation threshold. 
The discrete spikes are artifact of our treatment of relative veloc- 
ities, which are evaluated in a discrete way (Section [2.2t . If we 
consider a smooth distribution of the relative direction between the 
grains in collision, those spiky features should be smoothed out. 

To analyze the results in Fig. [2^, we examine the contribution 
from accretion and coagulation. We show the results for only accre- 
tion and only coagulation in Figs.|2j5 and c, respectively. Compar- 
ing all the panels in Fig. |2l we see that the evolution of grain size 
distribution is first driven by both accretion and coagulation, and is 
later caused by coagulation. In fact, the accretion time-scale is es- 
timated by equation ^ as a/a — r(a)/5(t), which is ~ 5 Myr 
for a ~ 0.001 [im if we use the value of ^ at t = (^ — 
0.3). Coagulation of small grains also occu rs on a similar time- 
scale. Indeed, according to equation (A5) in [Hirashita & Omukail 
(2009), the coagulation time-scale can be estimated as tcoag ^ 
4(nH/103 cm-3)-i(a/10-3 ^mf/^'is/S g cm-3)3/2(D/3 x 
10"^)"^(T'/10 K)"^/^ Myr, where V is the dust-to-gas ratio. 
Although coagulation is as effective as accretion in depleting the 
small-sized grains, it is the role of accretion that creates the strong 
bump around a ~ 0.002 \im. As we observe in Fig. I2J5, accretion 
stops at several Myr because of the depletion of gas phase met- 
als; thus, the evolution is purely driven by coagulation after several 
Myr. Comparing Figs.[2ji and c, the typical grain size achieved by 
coagulation is not sensitive to the presence of accretion, but the 
bump around ~ 0.02 \Lm is stronger in the presence of accretion. 
Thus, accretion helps to enhance the bump, while the typical grain 
size reached by coagulation is not affected by accretion. 

Since grain growth predominantly occurs at small sizes, it is 
expected that grain growth is sensitive to amin- In Fig. [3] we show 
the evolution of grain size distribution for Model B, in which we as- 
sume amin ~ 10 nm instead of 3 nm. In this model, grain growth by 
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Table 1. Adopted quantities. 



Species 


X 




mx [amu] ^ 


(X/H)0 


s [g cm ^] 


Silicate 


Si 


0.166 


28.1 


3.55 X 10-5 


3.3 


Graphite 


C 


1 


12 


3.63 X 10"* 


2.26 



^ For silicate, we assume a composition of Mgi.iFe().gSi04 iOraine & Led 
Il984) . 

° The atomic masses and the abundances are taken fromlCo>Jj200(ll) . 
The material densities are taken from lPraine & Le3 jl984h . 




Figure 1. Grain size distributions in Model A at t = 10 Myr. The thick solid, dashed, and dot-dashed lines show the results with both accretion and 
coagulation, only accretion, and only coagulation, respectively. The dotted line and the thin solid line present the initial size distribution and the analytic 
solution for accretion, respectively. Panels (a) and (b) show silicate and graphite, respectively. 



(a) Silioato (acc+coag) 




0.001 0.010 
a [/.Am] 



0.100 



10" 



a: 10" 



10" 



10" 



(b) Silicate (acc) 



initial 

1 Myr 
3 Myr 
10 Myr 
30 Myr 



0.001 0.010 
a [/.Am] 



0.100 



10" 



I 10" 



(c) Silicate (coag) 



10" 



10" 




0.001 



0.010 
a [/.im] 



0.100 



Figure 2. Same as Fig. [T] but at various t. The thick sohd, thick dashed, dot-dashed, thin solid, and thin dashed lines show the grain size distributions at 
t = 1, 3, 10, 30, and 100 Myr, respectively (in Panel b, since the grain size distributions are the same for t ^ 10 Myr, we do not show the size distribution at 
t = 100 Myr). Each panel shows the case with (a) both accretion and coagulation, (b) only accretion, and (c) only coagulation. 
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Figure 3. Same as Fig. but for a larger minimum grain radius, amin = 
1 rmi (Model B). The thick solid, thick dashed, dot-dashed, thin solid, and 
thin dashed lines show the cases for 4 = 1,3, 10, 30, and 100 Myr, respec- 
tively. 



accretion lasts 3 times longer than in Model A, because of 3 times 
larger amin (note that the grain growth time-scale is proportional 
to a; equation [5](- However, the growth of large grains (typically 
a > 0.003 \im) is not affected by the change of amin- 

Regardless of flmin, turbulence has a significant impact on 
the grain size distribution at a > 0.003 \im through coagulation, 
especially after accretion is terminated by th e depletion of gas- 
phase metals. In this paper, we have foll owed Orme l et al. 
for the grain motion driven by turbulence. lYari, Lazarian, & Dr ains 
( l2004h show that the grain motion in magnetized turbulence can 
also be excited by gyroresonance. However, the grain motion in 
dense medium is dominated by hydrodynamical drag even in their 
models, and the velocities that they calculated are similar to ours 
for a < a few x 10~^ (xm. 

In order to examine the significance of turbulent motion, we 
also calculate the case where the turbulent velocity component is 
neglected; that is, v — Uth (Model C). We show the result in Fig.|4] 
Comparing Fig.|4]with Fig.[2^, we observe that the effect of turbu- 
lence on the grain size distribution appears after f ~ 10 Myr, when 
coagulation begins to be important for a > 0.006 |J.m; for this grain 
size range, the grain motion is predominantly driven by turbulence. 
For t < IQ Myr, however, turbulence does not significantly affect 
the grain size distribution. 



3.2 Effect of coagulation on the grain mass growth 

Grain growth by accretion in interstellar clouds is one of the most 
important processes that govern the grain mass budget in an entire 
galactic system. Following HKl 1, we define the increased fraction 
of dust mass in a cloud, /3, as 



Figure 4. Same as Fig.|2ji but we show the case of no turbulence (Model C). 

where t^i is the lifetime of clouds hosting the grain growth[f| The 
dust mass in the cloud becomes (/? + 1) times as much as the initial 
value at the cloud lifetime, when the dust grown in the cloud returns 
in the diffuse ISM. By using /3, the increasing rate of dust mass in 
a galaxy should be written as 



dMd 



At 



(22) 



/3 = 



Pd(Tcl) 
Pd(0) 



(21) 



where i\/dust is the total dust mass in the galaxy, and X^i is the 
mass ratio of clouds hosting the grain growth to the total gas mass 
(HKll). 

Noting that the time-scales of accretion and coagulation are 
both scaled with the density of metals [oc (nu_Z)~^}, j3 can be 
regarded as a function of TcifinZ. In Fig. |5] we show j3 for Mod- 
els A and B to examine the dependence on grain size distribution 
(cimin = 0.3 nm and 1 nm, respectively). Because accretion occurs 
more efficiently at smaller grain sizes, /3 is larger for Omin = 0.3 
nm than for amin = 1 nm. Accretion saturates if the metals in gas 
phase are used up. Thus, /3 — >■ 5(0)/[l — C(0)] if the cloud lifetime 
is long enough. 

In Fig.[5]we also plot the cases where we neglect coagulation 
(dotted lines)[j There is a slight indication that coagulation sup- 
presses accretion. This is explained as follows: Coagulation pushes 
grains to larger sizes with the total volume conserved so that the 
surface-to-volume ratio of the grains decreases. Note that accre- 
tion rate is predominantly regulated by the surface-to-volume ra- 
tio. Nevertheless this suppression effect is very small and we can 
conclude that coagulation has little impact on accretion. This also 
means that we can neglect coagulation as long as we are interested 

In HKll, only accretion is considered, so is proportional to the mean 
value of . However, if we take coagulation into account, the grain number 
density also changes. Thus, in this paper, we cannot write /3 by using the 
mean value of . 

* We have confirmed that our results without coagulation matches the ana- 
lytical results by HKll. 
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Figure 5. Increment of dust mass, /3, defined in equation j2U as a function 
of the cloud lifetime t^i . To clarify the dependence on metallicity and den- 
sity, the scaling with Z and njj is also shown in the horizontal axis. The 
difference in the upper and lower lines with the same line species is the 
value of aj„in (0.3 nm for the upper lines and 1 nm for the lower fines). Tfie 
solid, dotted, and dasfied lines show the results with accretion and coagula- 
tion in the presence of turbulence, without coagulation (i.e. only accretion), 
and without turbulence (but with both accretion and coagulation consid- 
ered), respectively. Note that the upper solid line, the upper dashed line, 
the lower solid line, and the lower dashed fines are equivalent to Models 
A, C, B, and D, respectively. In the small window, we zoom in to clarify 
the distinction between the lines (vertical axis is on linear scale instead of 
logarithmic scale). 



in the total dust mass in a galaxy, although we should keep in mind 
that coagulation can affect the grain size distribution (Sec tion r3.il ). 

Finally, we also show the case where turbulence is neglected 
(Models C and D) in Fig.|5] We observe that the effect of turbulence 
on (5 is negligibly small. If flmin is larger, coagulation is very ineffi- 
cient if there is no turbulent motion; thus, the grain growth without 
turbulence is practically the same as the grain growth without co- 
agulation. We can also conclude that turbulence has little impact 
on the grain mass increase in clouds, although we should note that 
turbulence can affect the grain size distribution (Section FS.lt . 



3.3 Dependence on other parameters 

We briefly discuss the influence of various parameters. The depen- 
dence of the accretion time-scale on physical quantities are clar- 
ified in equations l |19l ) and MQ\ . Both accretion and coagulation 
follow scaling of time-scale oc {Znu)~^, so that we get the same 
grain size distribution at the same tZnu- Precisely, the grain veloc- 
ities driven by turbulence also depends on nn, but this dependence 
is only oc n^"^. The time-scale of grain growth is also inversely 
proportional to the sticking efficiency. Regarding the gas tempera- 
ture Tgas, both accretion and coagulation are again regulated by the 

— 1/2 

same scaling oc Tgas if coagulation is purely driven by thermal 
motion. Nevertheless, the gas temperature in cold gas is at most 



- 100 K jWilson. Walker. & Thomlevl ll997h : thus, the tempera- 
ture just causes a factor of 3 difference in the time-scales. Coag- 
ulation driven by turbulence is little affected by Tgas, because the 
dependence of the turbulence-driven grain velocity on gas temper- 
ature is weak (oc Tg/s'*). 

As mentioned in Section 12.31 we assume that r — 3.5. If 
r is smaller/larger, the number of grains at the smallest sizes is 
smaller/larger. Thus, as already shown in HKll (see their Fig. 3), 
accretion occurs more slowly/quickly. Because coagulation also oc- 
curs at the smallest grain sizes, the relative role of accretion and co- 
agulation does not change even if we change r. We should note that 
r = 3.5 is supported fr om the extinction curves in th e Milky Way 
and Magellanic Clouds jMathis et anil977l : |Peilll992h . Some theo- 
retical studies also indicate that disruption or shattering of grains in 
the diffuse ISM pr ocesses the gr ain size distribution to a power-law 
with r ~ 3.5 (e.g. lHellveilll97(lh . 

As shown above, grains at a '-^ fflmax are not affected by ac- 
cretion and coagulation. Thus, changing Omax does not affect our 
results. This is because both accretion and coagulation occurs effi- 
ciently for grains which have a large surface-to- volume ratio. 



4 DISCUSSION 

Now we discuss two issues related to accretion and coagulation. 
One is effects of these processes on the extinction curve, and the 
other is an implication for galaxy evolution. 



4.1 Effects on the extinction curve 



Following iHirashita & Ym\ ( l2009h and lO'Donnell & Mathi j 
(1997), we exa mine how the grain s ize evolution affects the 
extinction curve. IHirashita & "^ ( |2009|) adopt an initial condition 
which reproduces the Milky Way extinction curve, and incorporate 
the change of the grain size distribution by shattering and coag- 
ulation to examine how the extinction curve is modified by these 
processes. Here we examine effects of accretion and coagulation 
on the extinction curve. 

The initial condition is set up so that the grain abundanc e is 
consis tent with the mean extinction curve of the Milky Way bv lPeil 
(1992); that is, we adopt the same abundances of C and Si as listed 
in Table[T]but apply ^(0) — 0.25 and 0.15 for silicate and graphite, 
respectively, under Model A. This rough agreement is sufficient for 
our aim, since our aim is not detailed fitting of the extinction curve. 
We adopt the same dust optical propertie s and c alculation method 
of extinction curves as in IHirashita &YaiJj2009t) : The grain extinc- 
tion cross section as a function of grain size is derived from the Mie 
theory, and is weighted with the grain size distribution to obtain an 
extinction curve. 

In Fig. |6l we show the calculated extinction curves in units 
of magnitude per hydrogen. We present the results with both ac- 
cretion and coagulation, only accretion, and only coagulation at 
t = 10 Myr. First, we observe that accretion increases the grain 
opacity at all wavelengths because the grain mass grows. In or- 
der to examine the wavelength dependence, we also show the ra- 
tio between the extinction curve at t = 10 Myr to that at t = 
Myr (initial) in the lower panel. We find that accretion 'steepens' 
the extinction curve rather than 'flattens' it. Naively it may be ex- 
pected that the extinction curve would be flattened after accretion, 
because the mean grain size becomes large. Contrary to this expec- 
tation, the extinction curve becomes steeper after accretion. This is 
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1.0 0.5 



X [Mm] 
0.3 0.2 0.15 0.12 



Table 3. Summaiy of effects on the extinction curve. 
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Figure 6. Upper panel: Extinction curves in units of magnitude per hydro- 
gen. The dotted line (almost identical to the dot-dashed line) is the initial 
extinction curve before accretion and coagulation. The solid, dashed, and 
dot-dashed Unes show the results with both accretion and coagulation, only 
accretion, and only coagulation, respectively, at t = 10 Myr, based on the 
parameter sets in Mo del A. The filled squares represent the Galactic extinc- 
tion data taken from IPeil jl992h . Lower Panel: The extinction divided by 
the initial extinction. The correspondence between the line types and the 
models are the same as above. 



explained as follows. The extinction at a wavelength A is propor- 
tional to a"^ if a <C A, while it is proport ional to to a'^ (i.e., less sen - 
sitive to the grain size) if a > A/(27r) jSohren & Huffmarj|T983h . 
In other words, the ultraviolet (UV) extinction is more sensitive to 
the enhancement of small grains than the extinction at longer wave- 
lengths. Therefore, the extinction at shorter wavelengths increases 
more sensitively as a result of accretion. The 0.22- (j.m bump, which 
is due to small graphite, also becomes more prominent after accre- 
tion. At optical and near-infrared wavelengths, the extinction is rel- 
atively insensitive compared to other wavelengths, since it is more 
affected by the largest grains intact after accretion and coagulation. 
As A becomes large at mid-infrared and longer wavelengths (i.e. 
A ^ 271-ainax), the extinction is just proportional to the grain mass 
(or volume). 

As shown in Fig. |6] (comparison between the solid and dashed 
lines), coagulation flattens the UV extinction curve and lowers the 
carbon bump. This is because of the depletion of small grains. Co- 
agulation conserves the total grain mass (or volume); thus, the ex- 
tinction curve at A ^ 27ramax is unaffected by coagulation. 

Observations of the Milky Way show that extinction curves 
are flatter in the directions of dense clouds jMathis & Wallenhorsj 
ll98lHCardem et al.ll 19891) . This is not consistent with accretion, 
but is consistent with coagulation. ICardellietalJj 19891) also show 
that the extinction per hydrogen nucleus is smaller toward dense 
clouds than in the diffuse ISM. This requires strong coagulation. 
The cloud lifetime is possibly larger than 10 M yr, so that coagula- 
tion can push the grains to large sizes. Indeed, iKoda et al. I l l2009h 
suggest a lifetime of ^ 100 Myr or longer. Another possibility 
of strong coagulation is that coagulation additionally occurs in ex- 
tremely dense regions, e.g. regions associated with star formation 
jOrmel et alj2009l : lHirashita & Omukaill2009h . 



Process UV slope'' 


Carbon bump'' 


ref^ 


Accretion + 


+ 


1 


Coagulation — 




1, 2 


Shattering + 


+ 


2 


Shock destruction — 




3 



^ +: steeper; — : flatter 

^ +: stronger; — : weaker. 

1) this paper; 2) lHu-ashita & YanI {lOGfk : S'l lHirashita etall j200^ 



4.2 Significance in galaxy evolution 

Grain growth by accretion is shown to be important not only in 
nearby galaxies, but also in some high-redshift galaxy populations. 
As sh own in previous studies (HKll; llnouel I2OI2I ; lAsano et al.l 
l2012h . grain growth by accretion becomes prominent if the metal- 
licity exceeds a critical value. Those studies pointed out general 
importance of grain growth by accretion in nearby galaxies. It 
is also indicated that grain growth is indeed governing the dust 
abund ance in distant quasar s ^Michal owski et alj201Cll ; |Pipino et aF 



Gall et al 



201 1| ; IValianteetai]|201ll ; Asano et al.l l2012h (biit 

I2OII ). In these objects, it is expected that coagulation is also occur- 
ring, but according to our results in this paper, the effect of coagula- 
tion on the total dust content can be neglected. If we are interested 
in the evolution of grain size distribution, however, we should take 
coagulation into account, especially after a significant fraction of 
gas-phase metals are locked into dust grains. 

Finally, in Table [3] we briefly overview the effects of var- 
ious interstellar processes on the extinction curve, including 
those treated in our previous papers. The carbon bump and the 
UV slope are enhanced by shattering ( Hirashita & Yan 200^ 
and ac cretion (this paper), w hile they are flattened by coagu- 
lation (lHirashita& Yaj|2009l; this paper) and shock destruction 
jNozawa et alj|2007l; iHirashita et alj|2008h . We should consider at 
least these effects as processes determining the shape of the ex- 
tinction curves observed in a variety of gal axies including high- 
redshift quasars and gamma-ray bursts (e.g. iGallerani et al]|201(]| ; 
Ijang et al.ll201 ll , for recent observations). 



5 CONCLUSION 

In this paper, we have focused on grain growth processes in in- 
terstellar clouds. We have formulated and calculated the evolution 
of grain size distribution by accretion and coagulation in an in- 
terstellar cloud. Destructive processes are neglected in this paper, 
since they are considered to be ineffective in dense regions. We 
have confirmed the previous analytic results that grains smaller than 
~ 0.001 \im are depleted by accretion on a time-scale of several 
Myr. Coagulation also occurs on a similar time-scale, but accre- 
tion has a prominent impact in making a large bump in the grain 
size distribution around a ~ 0.001 \xm because of the grain mass 
increase. After a major part of the gas phase metals are used up, ac- 
cretion stops and coagulation is the only mechanism that pushes the 
grains toward larger sizes. It is only at this stage that grain motion 
driven by turbulence can affect the evolution of grain size distribu- 
tion. Regardless of whether accretion takes place or not, the typical 
grain size achieved by coagulation is similar. 

We have also examined the grain mass increase during the 
cloud lifetime by introducing /3 in equation | |2U such that the grain 
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mass becomes (1 + /?) times as much as the initial value. We have 
found that coagulation slightly suppresses accretion although this 
effect is not significant and can be neglected in discussing the evo- 
lution of the total grain mass in galaxies. In other words, coagula- 
tion can be neglected in discussing the evolution of the total grain 
mass in galaxies. 

Finally, we have investigated the effects of accretion and co- 
agulation on the extinction curve. We have found that accretion en- 
hances the UV slope and the carbon bump in spite of the increases 
in the mean grain radius. As expected, coagulation flattens these 
features. 
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APPENDIX A: DISCRETE FORMULATION 

Ror numerical calculation, we consider A*' discrete grain radii, and 
denote the lower and upper bounds of the ith (i = 1, ■ • ■ , A'^) bin 
as a'^l^ and a^^\ respectively. We adopt a'*^' = a'^^'jiS, Oq''' ~ '^i' 
and a^' — Uu with log 5 — (l/N) log(au/ai). We represent the 
grain radius and mass in the ith bin with Ui = {a^^-^ + a[^^ ) /2 
and nii = (47r/3)ai s. The boundary of the mass bin is defined as 
mf'^ = (47r/3)[o^''']'^s. The interval of logarithmic mass grids is 
denoted as = 35. We also discretize the time as t„ — nAt. 
We use integer indexes i and n to specify the discrete grain size 
and time, respectively. We adopt A'^ — 512, ai — 2 x 10"* \im, 
and Ou — 0.3 |J.m. We apply n(ao, t) — n{aN, t) = for the 
boundary condition. 
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hi Accretion 

Now we explain how to discretize equation l llOt . We denote the 
value of quantity Q at a discrete grid as Q" (recall that i and n 
specify grain-size and temporal grids, respectively). Then, we ob- 
tain the following equation as a discrete version of equation dlOt : 

^Ti + l n n -i 



At Am =3M.<^,. (A1) 

where the difference is evaluated based on upwind differencing. 
In hydrodynamical simulations, more elaborate differencing meth- 
ods have been developed, but the simple difference is enough 
for the grain size distribution (for example, we are not interested 
in 'shocks' or 'instabilities' when we are treating the grain size 
distribution in this paper). We apply Courant condition At < 
A/i/ maxi(/i") for every time step. From equation Jilt , we ob- 
serve that the minimum value of /i is realized at the maximum value 
of r(m), which is small when m is small (equation [Sjl. Therefore, 
the above Courant condition is practically evaluated at j = 1. 
Solving equation iAli for cr"+^, we obtain 

.„At At .A At „ 



A2 Coagulation 

The mass density of grains contained in the ith bin, p" (n is the 
time step), is defined as p" = a"miAfi. The time evolution of pi 
by coagulation can be written as 



--n+l ~n 
Pi -Pi 

At 



N N N 



i=i j=i 1=1 

(A3) 



and 



itVlk < Uc. 



aek = mkmt (A4) 
otherwise. 



where /3 is the sticking probability for coagulation, and the co- 
agulated mass m*jag(j) is determined as follows: m^ja^g(i) = 
m£ if Jn'^j ^ me + nij < m''^'[f| otherwise rr/Joa.g.ii) = 
0. Coagulation is assumed to occur only if the relative veloc- 
ity is less than the coagulati on thr e shold velocity t'co ae- based 
on IChokshi. Tielens . & Holle nbaclj h993h. j Pominik & Tielen j 



( Il997l) . and lYanetal.1 ll2004) (see 'Hirashita & Yani2009l , for fur- 



ther details). The cross-section for the coagulation is aik = n{ai -\- 
o-k)^. We assume /3 — 1 for the sticking probability. 

To confirm the v alidity of our nu merical scheme of coagu- 
lation, we checked in iHirashita & ^^[ 2009 ) that the total dust 
mass is conserved. We also observe in Hirashita & Omukail j2009h 
that coagulation proceeds on an analytically estimated time-scale: 
this check is roughly equivalent with that in the appendix D of 
iBrauer, DuUemond. & Hennind feOOSi) . 
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